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Abstract 

A system of two interacting photon modes, without constraints on the photon number, in the 
presence of a Kerr nonlinearity, exhibits BEC if the transfer amplitude is greater than the mode 
frequency. A symmetry-breaking field (SBF) can be introduced by taking into account a classical 
electron current. The ground state, in the limit of small nonlinearity, becomes a squeezed state, 
and thus the modes become entangled. The smaller is the SBF, the greater is entanglement. 
Superfluid-like behavior is observed in the study of entanglement growth from an initial coherent 
state, since in the short-time range the growth does not depend on the SBF amplitude, and on 
the initial state amplitude. On the other hand, the latter is the only parameter which determines 
entanglement in the absence of the SBF. 

PACS numbers: 03.75.Gg 
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Quantum entanglement is considered as a fundamental resource in quantum information 
science In the last years, quantum information with continuous variables (CVs) 

aa 

has receveid a lot of attention because of the simplicity of preparing, unitarily manipulating, 
and measuring quantum states. Various CV applications have been considered such as 
teleportation [4j and cryptography [5(. 

At the same time, quantum phase transitions (QPTs) are recognized as a fundamental 
concept in quantum many-body systems {(J. In particular, the transition to a ground state 
with a macroscopic occupation, the Bose-Einstein condensation (BEC) regime, is widely 
studied both theoretically and experimentally Recently, a great effort has been 

devoted in understanding how quantum correlations (such as entanglement) are related to 



the critical behavior of systems exhibiting QPT effects [10, 
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Remarkable examples of QPTs in CV systems are superradiance 
Einstein condensation (BEC) in microcavity polaritons [14J ]. Some years ago, Chiao in- 
vestigated the possibility of observing BEC and superfluidity in photon gases [la ] . In fact, 
since photon quantized in a three-dimension cavity are massless and their chemical potential 
is zero, BEC transition would appear forbidden. To overcome this difficulty, Chiao proposed 
a two-dimensional array of Fabry-Perot cavities where an effective chemical potential can 
be defined due to quantization conditions. A similar approach has been also proposed by 
Navez fl- 

In this paper we show that a ground state with a finite occupation (BEC) for pure pho- 
tonic systems is indeed possible without the need of special architectures. We will consider 
two photon modes, coupled via an exchange interaction, in the presence of a weak nonlinear- 
ity The emergence of an effective "chemical potential" , manifests itself as a consequence of 
the unitary transformation which diagonalizes the quadratic part of the Hamiltonian. The 
normal phase corresponds to the situation where the photon frequency exceeds the interac- 
tion coupling, while in the opposite case a BEC phase comes out. In the anomalous phase, 
a finite amount of entanglement in the ground state is predicted between the two modes 
which grows with the occupation number. In Ref. [l3] a two- mode boson entanglement 
across the transition from a normal to a superradiant phase is discussed. Squeezing and 
entanglement in the ground state are present in both phases. The main difference with our 
model is that squeezing and entanglement appear only in the anomalous phase. Other sev- 
eral model, discussing the interaction of two e.m. modes in nonlinear media have been also 
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considered [l^, Q- While ground-state properties allow the investigation of the quantum 
nature of the system, quantum optics experiments are usually set to observe the transient 
regime determined by cavity losses Q. Then, it seems important the discussion of how an 
initially factorized state can develop entanglement in its time evolution. The evolution of 
entanglement from a given initial state has been studied in Ref. [21| for Fock and coherent 
states. In the cases considered there, the evolution does not depend on the presence of BEC. 
On the other hand, if a symmetry-breaking field (SBF), physically achievable through the 
interaction of the e.m. radiation with a classical electron current is present, a "superfluid- 
like" phase appears. The main difference with the work of Ref. 2l| is the existence of an 
initial time range where the entanglement growth does not depend on the SBF amplitude. 

Let us start with the description of the model by considering two coupled photon systems 
dj (i = 1, 2) with the same energy frequency uj, a transfer interaction of strength w, and a 
nonlinear Kerr interaction of strength g 

H = oj (^a\ai + a\a-2^ — w (a\a 2 + a\a^j + g (ri\ + n 2 ) 2 (1) 

The two systems could be associated to the polarization components of the same e.m. field. 
This is the case where the assumption about the interaction is sensible. 

To diagonalize H we perform a rotation on the original degrees of freedom which defines 
new variables a — (a% + a 2 ) / \p2 and @ — (a\ — a 2 ) / \/2. 

This unitary transformation cancels the transfer term 

H = (u - w) a ] a + (u + w) ftp + g (a ] a + ft[3) 2 . (2) 

We observe that the last term, even if small, becomes important for w > u, when the 
vacuum photon state is no more the ground state. In fact, in the absence of interaction, the 
Hamiltonian is no more bounded from below. 

The Hamiltonian (J2]) is diagonal in the basis spanned by the Fock states in the (at, (3) 



2l|. The ground state is the vacuum when (u — w) > 0, whereas it becomes 
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\n a ) \0p) when (uj — w) < , where n a is the integer closest to (w — uj) / (2g). In the 
latter case the corresponding eigenvalue is (u — w)n a + gn 2 a . Ground state properties are 
remarkably different in the two cases. Studying the entanglement of formation between the 
modes a and b by means of the Von Neumann entropy S of either of the reduced states, 
two regimes appear. In the normal phase, when the ground state is the vacuum, there is no 



entanglement, while when the ground state is \n a ) \0p) the following expression for S holds 
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which shows how the Von Neumann entropy grows with n a . In Fig. [T]the behavior of S(a, b) 
is plotted as a function of the critical parameter w/(u + g). 
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FIG. 1: Von Neumann entropy of the subsystems a and b versus the critical parameter w/(u) + g). 
In the normal phase no entanglement is detected, while beyond the critical point the ground state 
is no longer factorizable. The system's parameter are u = 1 and g = 0.01. 

As said, ground state properties are significant, but more important can be the achieve- 
ment of an experimental method to create entanglement and to exploit it in quantum in- 



formation protocols. The authors of Ref. 2l(] have studied the entanglement generation 



induced in the time evolution of Fock states and coherent state. Let us focalize our atten- 
tion on these latter states. We note that the dynamical entanglement does not depend on 
the presence or the absence of the BEC phase. This result can be understood by observing 
that uj can be dropped through the introduction of a "local" unitary transformation which 
cannot determine any change in the entanglement value. Furthermore, the amplitude of 
the incoming coherent state influences deeply the degree of entanglement. Then, low in- 
tensity two-mode laser fields are expected to develop a little amount of entanglement. The 
introduction of a finite SBF modifies these properties in a significant way. Let consider 
a field of amplitude X/V2 acting on both the modes a and b . The new Hamiltonian is 
H\ = H — A (cJ + a) . This field does not represent a purely mathematical tool introduced 



to describe the emergence of a superfTuid phase, as usual in boson particles systems, but 
it can be physically realized through the interaction of the e.m. modes with a "classical" 



non- fluctuating electron current [20(], and, then, it can be kept finite. An exact solution for 
this model does not exist anymore, and a perturbation theory is worth to be developed. 
First, in order to study the properties of the ground state of H, it is easy to recognize 
that the vacuum state for the mode f3 minimizes the energy. Then, we write the ground 
state \G) as the product \G a ) \0p)- Let us try to do a perturbation theory starting from the 
solution known for A = 0. In the normal phase the correction to the ground state energy 
is AE = —A 2 / (u — w + g). The result is substantially different in the anomalous phase. 
When (w — u) /2g is an integer, the first excited states are degenerate (\n a + 1) and \n a — 1) 
have the same energy) and the energy gap is g: AE = —A 2 (2n a + 1) / g. If (w — uj) /2g is 
not integer, the solution is just slightly different. A qualitatively different solution arise 
when (w — u) /2g falls exactly between two integer numbers. Now the ground state is de- 
generate, and energy corrections are of order A. Apart from these special values, in the limit 
of small nonlinearity, which corresponds to g < A, the energy spreading is greater than the 
gap. Then, the system becomes unstable with respect to this perturbation, and a different 
approach is in order. 

Next, we show that the effect of nonlinearity can be taken into account introducing a 
"superfluid-like" phase in the limit of small g. We mean that in the new ground state the 
annihilation operator a has a nonvanishing mean value v. It is then natural to introduce a 
translation of amplitude v (a — > a + v). The Hamiltonian becomes H\ = Eq + Hi + Hi + 
H 3 + if 4 , where 

E = (u - w + g) v 2 + gv 4 - 2Xu, (4) 
Hi = (a j + a) [(u-w + g)u + 2gu 3 - A] , (5) 
H 2 = (lu - w + g) n a + (u + w + g) rip + gv 2 [An a + a t2 + a 2 ] , (6) 
i7 3 = 2gu [(cr + a) np + a^n a + n a a] , (7) 
H 4 = g [a ]2 a 2 + fi 2 (3 2 + 2n p n a ] . (8) 

The weak-coupling approximation amounts to take into account terms of order gv 2 , and 
to neglect terms of order yjg and g. In this limit we can disregard H 3 and H4. The 
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condensate amplitude v can be fixed by minimizing Eq\ 

(Lu-w + g)v + 2gu 3 - A = 0. (9) 

We note that this conditions makes Hi vanishing. 

Two different behaviors appear depending on the sign of the coefficient of v. If 
(a? — w + g) > we have v ~ A/ (uo — w + g). In this case we have a "normal" behav- 
ior, in the sense that condensate amplitude vanishes with A. On the other hand, for 
(u — w + g) < 0, a solution for the minimum of E (v) is given by v ~ v* + A/ (4gz/* 2 ), 



with i/* = \j w 2 9 9 ■ This solution corresponds to a superfluid state. 

Taking into account Eq.fjSJ), we obtain for the Hamiltonian if 2 = H a + Hp, where 

H a = -n a + gv 2 (2n a + a ]2 + a 2 ) (10) 
if/3 = (u + w + g) rip (11) 

Qualitatively, we observe that in the normal phase H a is dominated by the term (\/v)n a 
and the ground state is expected to be very close to the vacuum state. Conversely, in 
the anomalous phase a squeezed state on the mode a will appear. Thus, entanglement is 
expected to come out between the modes a and b, in analogy with the well known effect 
concerning the passage of squeezed light through a beam-splitter. 

It is worth noting that for A = a quadratic dependence on the operator x a = 
(oft + a) I \p2 comes out. Thus 0C> q IS 3j constant of motion, and the canonical conjugate 
variable p a = i (oft — eft) / \/2 freely diffuses. This feature has been recognized in Bose par- 
ticles systems by Lewenstein 22]. It is commonly assumed (see for a review 23j) that 



corrections of order g are able to remove this singular behavior. This result is known in the 



241 ]. and rejected on 



Bose particles system as the Hartree-Fock-Bogoliubov approximation 
the basis of general arguments (Goldstone theorem). 

For finite, even if small A, it exists a canonical transformation (the Bogoliubov transfor- 
mation) which diagonalizes H a . If we define 7 = cosh^a — sinh^a"'', we get H a = e^j + eo, 



where e = y^(^ + 4g^ 2 ) and eo is the zero-point energy. The parameter 9 is defined 
by tanh26 l = —2gu 2 / (^ + 2gu 2 ). The ground state \G a ) of H a is the squeezed state 
\G a ) = exp [| (aft aft — otaftft |0). In the original space a, b, taking into account of the trans- 
lation performed on a, 

'6 



\G a ,b) = exp 



- ( aft aft + b ] b ] + 2aftb ] - \[2v (a f + 6^)^ — h.c. 



|0,0>. (12) 
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In the Schrodinger representation, where the coordinates are x a = (al + a) / V% and Xb 
(&t + b) I y/2, the wave function reads as 



(x a , x b ) = exp 



-2B 



1 



e 2 

[{Xa + X b ) + 2u] - ~{X a 



Xb 



(13) 



Given the wave function, and then a complete description of the state, we characterize the 
phase transition through the entanglement of formation between the subsystems a, b. Since 
we have a pure two-mode state, the Von Neumann entropy of either of the reduced states is 
a unique entanglement measure. From our wave function, the Von Neumann entropy is 25] 



Of 0\ Of 0\ 

S (a, b) = cosh 2 - log I cosh 2 - j — sinh 2 - log I sinh 2 - J . (14) 
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FIG. 2: Von Neumann entropy of the subsystems a and b versus the critical parameter w/(uj + g) 
in the presence of a SBF A. The system's parameter are uj = 1 and g = 0.01. The orange (light 
gray) line refers to A = 0.3, while the black (dark gray) line corresponds to A = 0.1 

In Fig. (a, b) is plotted as a function of the critical parameter wj {uj + g) for different 
values of A. In the normal phase there is not entanglement, while in the condensate phase 
S (a, b) becomes drastically different from zero. The plot shows that the smaller is A the 
bigger is the amount of entanglement in the ground state. 

Coming back to the problem represented by a dynamical generation of entanglement 
starting by a factorized state, we next consider the dynamics of a coherent state |$) = |z/, v') 
with the same amplitude v' on both the modes a and b. Furthermore, we assume that the 
system is in the condensate phase. In the space of a and (3 this state corresponds to the 
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coherent state of amplitude z/a/2 for the mode a times the vacuum for (3. To solve the 
problem of calculating e~ lHt |<&) we use the trivial identities a (—t) |$ (t)) = u'y/2 |$ (i)) and 
/3 (— £) |$ (£)) = 0. The evolution of the annihilation operators in the Heisenberg picture 
can be obtained noting that 7 (t) = exp (— iet) 7 and (3 (t) = exp [— i (u + w + g) t] (3. After 
the substitution of these expressions and after some algebraic manipulations, the following 
result is established: 

a(t) = f(t)a + f (t)ai + h(t), (15) 

where 



/ (t) = cos et — i cosh 29 sin et, 
f {t) = i sinh 29 sin et, 



h{t) = v — v (cos et — ie 



-2B 



sm 



et) 



(16) 
(17) 
(18) 



Through these results, two coupled partial differential equations can be written for the 
wave function ty(x a ,Xb) = (%a, Xb\& (t)) which allow to solve the problem. The equations 
are 



(/ + /') (Xa + X b ) + (/- /') 



d d 
+ 



dx a dxb 



+ V2{h 



x b + 



_d_ 

OXn 



_d_ 

dx b 



* (x a ,x b , -t) 
* (x a ,x b ,t) 








(19) 
(20) 



Apart from a normalization constant, the solution is 



* (x a , x b , -t) = exp 



h-u' 



x exp 



/-/ 
-5 



- (x + x 6 ) 



(21) 



J-f'J Ka °' \f-f, 
The last term is responsible for entanglement, and is intrinsically related to the squeezing 
effect. In fact, it would be zero for 9 = 0. The initial amplitude v' does not play any role 
in the entanglement evolution. The Von Neumann entropy (whose analytical expression is 
rather intricate) is plotted in Fig. [3] as a function of time. 

An oscillation with frequency 2e is observed. As a crucial point of this derivation, we 
remark that the entanglement is completely independent from the amplitude of the incoming 
coherent state. This result is a specific signature of the superfluid phase. In fact, in the 
absence of the SBF, the evolution of a coherent state involves an amount of entanglement 
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FIG. 3: Dynamical evolution of entanglement as a function of time. The initial state is the product 
of two coherent states on the modes a and b. The system's parameter are uj = 1, w = 2, g = 0.1, 
and A = 0.11. 

which is strongly related to its amplitude and goes to zero as v' does. Further, the short- 
time evolution shows a weak dependence on the SBF amplitude. The limit of validity of the 
free- diffusion regime is given by A 2gu 3 together with et <C 1. 

In conclusion, we described the instability which occurs in a two-mode photon system in 
the regime where the internal coupling is stronger than the photon frequency. The ground 
state of the system is the vacuum below the critical point, while in the other region a 
finite population of photons appears. In the presence of a classical electron current, a 
"superffuid-like" phase manifests itself and gives rise to some interesting effects which could 
experimentally observed, such as the generation of an amount of entanglement in a time- 
evolution of factorized coherent states independently from the incoming radiation amplitude. 
Even if a zero- dimensional system a superffuid phase cannot be reached, by keeping the SBF 
to a finite value, we predict a short-time superffuid-like behavior. We thank S. Paganelli for 
continuous and very useful discussions. S. Solimeno, P. Mataloni, and F. Illuminati are also 
acknowledged. 
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